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Abstract 

In this paper, we consider the problem of Mutually Unbiased Bases in 
prime dimension d. It is known to provide exactly d + 1 mutually unbiased 
bases. We revisit this problem using a class of circulant d x d matrices. The 
constructive proof of a set of d + 1 mutually unbiased bases follows, together 
with a set of properties of Gauss sums, and of bi-unimodular sequences. 



1 INTRODUCTION 

Mutually Unbiased Bases (MUB's) are a set {Bq, Bn} of orthonormal bases of 
such that the scalar product in of any vector in Bj with any vector in 
^k, Vj 7^ A; is of modulus Starting from the natural base Bq consiting of 

vectors Vi = (1, 0, ...0), V2 = (0, 1, 0, 0), ...,Vd = (0, 0, 1), it is known that this 
problem reduces to find N unitary Hadamard matrices Pj such that PjPk is also a 
unitary Hadamard matrix Vj ^ k. (A unitary matrix is Hadamard if all its entries 
are of modulus d'^^"^). The problem has been solved in prime power dimension 
d = p", p,n G N,p prime, and yields exactly d + 1 MUB's which is the maximum 
number of MUB's ([2] and references herein contained). 

If d is factorizable in p^^p^^--- with pi ^ pj prime numbers, it is known that one 
has at least N = min p™' [5j. 

In this paper we show that in prime dimension d = p, the Discrete Fourier Transform 
F together with a suitable circulant matrix C allow to construct a set of d+1 MUB's. 
In addition this construction allows to obtain, as a by-product, a set of properties 
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of Gauss sums of the following form : 

= v^, Vj G ¥(i, I G coprime with d, d> 3 

(1.1) 

(Fd is the field of residues modulo d). A direct proof of this property is given in [12]. 
Similar results on generalized Gauss Sums appear in |1]. The definition of F and 
of circulant matrices is given below. The natural role played by circulant matrices 
in this context is a new result. The circulant unitary matrices are known to be in 
one-to-one correspondence with the bi-unimodular sequences c = (ci, C2, q) 
namely sequences such that \cj\ = \{Fc)j\ = 1, where F is the discrete Fourier trans- 
form. Not surprisingly Gauss sums appear naturally in this context, since suitable 
Gauss sequences are examples of bi-unimodular sequences. 

At the end of this paper, we consider the case of non-prime dimension, and show 
that Gauss sums properties can be deduced in the odd and in the even dimension 
cases. In a forthcoming work we shall consider the case of prime power dimensions 
and show that the theory of block-circulant matrices with circulant blocks solve the 
MUB problem in that case. 

A d X d matrix is Hadamard if all its entries have equal moduli [7] , and 

H*H = dl 



^ f2iTT \lk(k + l) 

k=0 ^ L -I ^ 



Definition 1.1 A dx d matrix H is Hadamard if \Hj^k\ is constant Vj, k = 1, ...,d, 
and H*H = d 1. We call H a "unitary Hadamard matrix" if 

d 

\Hj^k\ = d'^^^, yj,k = l,...,d, and ^ HjjHi^k = 5j,k 

1=1 

Definition 1.2 A dxd matrix C is ca/Zed circulant ^6^, and denoted ciTc{co, ci, Cd-i), 
if all its rows and columns are successive circular permutations of the first. It is of 
the form 





f Co 


Cd-1 ■ 


. cA 


c = 


Cl 


Co • 


• C2 




\Cd-l 


Cd-2 ■ 


• Co/ 



Proposition 1.3 (i) The set C of all d x d circulant matrices is a commutative 
algebra. 

(ii) C is a subset o/ normal matrices 

(ii) Let V = circ(0, 0, 1). Clearly = 1. Then C is circular if and only if 
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it commutes with V, and one has for any sequence c = (ci,...,Cd) G C^, C = 
circ(ci, ...,Crf) : 

C = P,{y) = Col + Cd-iV + ... + ciV"-^ 
where Pc is the polynomial 

d 
k=0 

Proof : See [6] 

For the product of two circulant matrices C, C (that therefore commute with V), 
one has 

VCC = CVC = CC'V 
which estabhshes that CC is indeed circulant. 

Moreover it is well-known that there is a close link between the circulant matrices 
and the discrete Fourier transform. Namely the latter diagonalizes all the circulant 
matrices. The Discrete Fourier Transform is defined by the following d x d unitary 
matrix F with matrix elements : 

F,,, = rf-V2exp(^), j,k = 0,l,...,d-l (1.2) 

Proposition 1.4 (i) The circulant matrix V = circ(0,0, 1) is such that 

F*VF =U = diag(l, uj, cu^"^) 

where 

^ = exp ^— J (1.3) 
(ii) Let C = circ(co, Ci, Cd-i) be a circtdant d x d matrix. Then 

F*CF = \/d diag(co, c_i, c_(rf_i)) 

where 



k=0 



Proof : (i) It is enough to check that Vk, the k-th column vector of F which has 
components 



(^^fc)i = j,k = 0,l,...,d-l 

is eigenvector of V with eigenvalue u^, which is immediate since 
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(ii) One has (Proposition 1.3(ii)) 



d-1 



C = J2ckV-' 

k=0 

Thus 

d-1 d-1 

F*CF=Y,Ck{F*VF)-'' = J2ckU-' = diagK,...,cid-i) 

fc=0 k=0 

with 

d-1 



k=0 

Lemma 1.5 For any k ^ N, we denote by [k] the rest of the division of k by d. 
Given any sequence c = (cq, G its autocorrelation function obeys 

d-1 d-1 
A:=0 1=0 

where the Fourier transform c of c has been defined in 
Proof : For any j = 0, ...,d — 1, one has 

d-1 ^ d-1 d-1 d-1 



k'l 

d 

1=0 1=0 k=0 k'=0 



since 



d-1 ^ d-1 d-1 

k,k'=0 1=0 k=0 

d-'J2^^^''~'~'^ = ^>^',iHk] 
1=0 

It is known ([1], [13]) that circulant unitary Hadamard matrices are in one-to- 
one correspondance with bi-unimodular sequences c = (cq, Ci, Cd-i)- 

Definition 1.6 A sequence c = (cq, Ci, is called bi-unimodular if one has 
\cj\ = \cj\ = 1, Vj = 1, ...,d, where Cj is defined by 



Proposition 1.7 Let (cq, Q-i) be a bi-unimodular sequence. Then the circulant 
matrix C = d~^^'^ciTc{co, ci, is an unitary Hadamard matrix. 
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Proof : This is a standard "if and only if" statement. One uses Lemma 1.5 : 

d-1 d-1 
k=0 1=0 

But since |q| = 1, Wl = 0, ...,d — 1, the RHS is simply d6jfi, and therefore 

d-1 

CkC[k+j] = d 6jfl 

k=0 

which proves the unitarity of the circulant Hadamard matrix C. 

In all that follows we call indifferently F or Pq the discrete Fourier transform. 

In [H], the authors introduce for any dimension d being the power of a prime number 
a set of operators called "rotation operators" which can be viewed as "circulant 
matrices" (this property is however not put forward explicitely by the authors). In 
this paper, restricting ourselves to the prime number case, we show that these 
operators can be used to define a set of d+1 Mutually Unbiased Bases in dimension d. 

Mutually Unbiased bases are extensively studied in the framework of Quantum 
Information Theory. They are defined as follows : 

Definition 1.8 A set {Bi,B2, ...,Bm} of orthonormal bases ofC^ is called MUB if 
for any vector h^^ G B^ and any h^^ ■* G Bk' one has 

\hf .hf\ = d'^l\ \/k^k' = l,...,m, Vj,/ = l,...,d 
where the dot represents the Hermitian scalar product in C^. 

Remark 1.9 It is trivial to show that if the orthonormal bases Bk are the column 
vectors of an unitary matrix A^, then the property that must satisfy the A^ 's in order 
that {1, Bi, Bfn} are MUB's is that all A^, k = l,...,m and AlA^', k ^ k' = 
l,...,m are unitary Hadamard matrices. Namely if uj, Vk are column vectors for 
unitary matrices A, A' respectively, then 

Thus unitary Hadamard matrices play a major role in the MUB problem. 

It is known that the maximum number of MUB's in any dimension d is d + 1, 
and that this number is attained if ci = p"*, p being a prime number. In this paper, 
restricting ourselves to m = 1, we revisit the proof of this property, using circulant 
matrices introduced by [8|. We then show that it implies beautiful properties of 
Gauss sums, namely the following ([12]) : 
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Proposition 1.10 Let d > 3 be an odd number. Then \/k 
with d the sequences 
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(k) 



exp 



iTckj{j + 1] 
d 



1, ...,d — 1 coprime 
(1.5) 



j=0,...,d-l 



are bi-unimodular. Thus U.l\) holds true. 



This property will appear as a subproduct of our study of MUB's for d a prime 
number > 3 via the circulant matrices method. As stressed above, the link between 
circulant matrices and bi-unimodular sequences is well established. What is new 
here is the fact that the MUB problem via a circulant matrix method allows to 
recover the bi-unimodularity of Gauss sequences. The crucial role played by the 
Gauss sequence is due to the crucial role played by the Discrete Fourier Transform 
(or in other therms the Fourier- Vandermonde matrices) in the MUB problem for 
prime numbers. Let us introduce it now explicitely. 

It is known since Schwinger ([H]) that a simple toolbox of unitary d x d matrices 
sometimes refered to as "generalized Pauli matrices" U, V allows to find MUB's. 
U, V generate the discrete Weyl-Heisenberg group [H]. Denote by u the primitive 
root of unity (II. 3p . The matrix U is simply 



U = diagfl, UJ, UJ 



which generalizes the Pauli matrix to dimensions higher than two. The matrix 
V generalizes ax ■ 



V = circ(0,0,...,l) 



/o 


1 


. 


• ^\ 








1 . 


. 








. 


. 1 


V 





. 





Then one has the following result : 

Theorem 1.11 (i) The U, V matrices obey the uj- commutation rule : 

VU = ujUV 

(a) The Discrete Fourier Transform matrix Pq = F defined by U.^) . namely 





(I 


1 


1 


1 \ 




1 


UJ 






"^'-Td 


1 


u^ 


U^ . 














V 
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diagonalizes V, namely 



(in) One has Pq = 1 



V = PqUPq = PoU*Po 



(ii) is simply Proposition 1.4(i). For the proof of (iii) reminiscent to the properties 
of continuous Fourier transform, it is enough to check that 



/I 




p2 



W 



10 
\0 1 



1 





0/ 



Thus W = lm dimension d = 2, and = 1, Vc/ > 3. 

In [5] we have shown that for d odd one can add to the general toolbox of unitary 
Schwinger matrices U, V a diagonal matrix D of the form 



D = diag(l, uj, 



k{k+l)/2 



1) 



SO that the MUB problem for odd prime dimension reduces to properties of U, V, D, 
and that certain properties of quadratic Gauss sums follow as a by-product. 



2 THE d=2 CASE 

We have f/ = cr^ and V = ax, cTz, o'x being the usual Pauli matrices. Since UV = ay, 
finding MUB's in dimension d = 2 amounts to diagonalize a^, ay. One has : 

a^ = PoCTzPq 
ay = Pia^Pi 



with 



Pi is circulant. 

Proposition 2.1 The set {l,Po,Pi} defines three MUB's in dimension 2. 

Since the matrices Pq, P\ are trivially unitary Hadamard matrices, it is enough 
to check that PqP\ is itself a unitary Hadamard matrix, which holds true since 



—i —i 



V2 
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3 THE PRIME DIMENSION d>3 

d being prime, we denote by the Galois field of integers mod d. 

Let us recall the definition of the "rotation operator" of [8], which, as already stressed 

in nothing but a circulant matrix in the odd prime dimension d. 

Definition 3.1 Define R as an unitary operator commuting with V and diagonal- 
izing VU. 

Proposition 3.2 (i) R is a circulant matrix, 
(a) R'' is also circulant V/c G Z. 

This follows from Proposition 1.3. 



Therefore we are led to consider a subclass of circulant matrices that are unitary. 

They must satisfy : 

Vfc = 0, d — 1 jcfcl = d'^^"^ and 

Vfc = 1, d — 1 Yl'jZo ^j^ld-k+j] = (orthogonality condition). 



Now it remains to show that such a matrix R exists. In [5] we have constructed 
a unitary matrix Pi that diagonalizes VU. It is defined as 

Pi = D-'Po 

for any d odd integer (not necessarily prime). 
We have established the following result : 

Proposition 3.3 (i) For any odd integer d, the matrix P^Pi is a unitary Hadamard 
matrix. 

(a) For d > 3 odd integer, let P^ '■= D^^Pq. Then P^Pk is a unitary Hadamard 

matrix for all k coprime with d. 

(Hi) The matrix D defined above is such that 

coprime with d 

For the simple proof of this result see [5]. (iii) is a simple consequence of (i) 
and(ii). Namely the eigenvector of VU belonging to the eigenvalue 1 has components 

Vj = d ' LU 2 

Since P^Pi is unitary Hadamard matrix, the element of PqPi of the first row and 
first column is simply 



3=0 



d-l 

UJ 2 
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and since its modulus must be d we obtain (iii) for k = 1. The proof for any k 
coprime with d can be obtained similarly using (ii). 

The "problem" is that Pi is not circulant. However the circulant matrix R is ob- 
tained from Pi by multiplying the kth column vector of Pi by a phase which is 

fc(fc-i) 

This operation preserves the fact that it is unitary and that it diagonalizes VU. We 
thus have : 

VU = PiUPi = RUR* (3.6) 
Proposition 3.4 Let R be the matrix : 

R=d-'/'ciTc (l,cu-^^-^...,^-^('=+^)/^...,l) 

It is a unitary Hadamard matrix. 

Proof : By construction it is a unitary Hadamard matrix, since Pi is. To prove 
the fact that it is circulant, it is enough to know that 

(Pi),,, = d-^'v-^ 

thus 

r> 7-1/2 ,7; 3(3+1) fc(fc-l) , (j-k)(j-k+l) 

sincG 

+ k(k-i) _ (j-k)(j-k+i) 

Thus all column vectors are obtained from the first by the circularity property. 
Furthermore the P*^ have the property that they diagonalize V'^U, V/c = 1, d—1: 

Theorem 3.5 (i) R = aP^DP* where a := d-^/'^Y^'J-^uj-^^^+^^l'^ is a complex 
number of modulus 1. 

(ii) R'^ = a'^l where 1 denotes the unity d x d matrix, 
(iii) 

R^U{R*f ^V^U, VA; = 0,...,d 

The proof is extremely simple : 

(i) We have shown that c, = uj 2 is a bi-unimodular sequence, thus a is a 
complex number of modulus one. Moreover from Proposition 1.4 (ii), the unitary 
matrix P = d~-^/^circ(cj) is such that 

P = Po*PPo = diag(c_fc) 
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But 

d-1 



1 \ ^ ih k{k+i) 



since 

d-l d-1 d-1 



j=0 j=0 i=0 

Therefore 



fc(fc+i) 

R = adiag(u; 2 ) = aD 

(ii) is simply a consequence of (i) since D'^ = 1. 

(iii) is obtained by recurrence. Namely it is true for k = 1 hj fl3.6p . For k > 2 one 
has : 

R^U{R*f = RV'^-^UR* = V'^'^RUR* 
since R commutes with V . But 

V^-^RUR* = V''-\VU) = V^U 

There is a direct link between the matrices that diagonalize VU^ and the R'^ 
that diagonalize V''U : 

Theorem 3.6 a being the complex number of modulus one defined above, we have 
for any k = 0, ...,d — 1 

Pu = a''P*R-''P^ 

Proof : In [5] we have proven that Pk = D~^Pq diagonalizes VU^. But one has 



j^-k ^ jjd-k ^ a^'dp.^d-kp^ 
so the result follows immediately. 

Corollary 3.7 For any k = 1, d—l,R^isa unitary Hadamard circulant matrix 
when d > 3 is prime. 

Proof : i?^ is circulant and unitary since R is. Therefore we have only to check 
that it is Hadamard. We have : 

R-' = a-'P^P*PkP^ 

But we recall that Pq equals the permutation matrix W. Thus all matrix elements 
of R'^ equal, up to a phase, some matrix elements of P^Pk- But we have established 
in [5j that the matrix PQPk is unitary Hadamard Vfc = l,...,(i — 1, thus all its matrix 
elements are equal in modulus to d~^/'^. This completes the proof. 



Now we show how this property of the matrix R reflects itself in Gauss sums 
properties. 
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Proposition 3.8 (i) Let d he an odd prime. Then for any k = 1,2, ...,d — 1 R'^ is 

an unitary Hadamard matrix if and only if one has 



^exp ( — [kf + j{k + 2m)] 

j=0 



Vd, Vm = —d + 1, ...,d — 1 



(a) Under the same conditions is a unitary Hadamard matrix if and only if 

k-l 



5Z ( T^-^^ + {k + 2m)j 

j=0 



Vk 



Proof : Since R'^ = PqD^ Pq , the matrix elements of R'' are 



k '^-1 



i=0 



Since R^ is circulant, it is unitary Hadamard matrix if and only if the matrix ele- 
ments of the first column (1=0) are of modulus d~^^'^, thus if and only if 



^exp (— [A;j2+j(A; + 2m)] 
j=o 



which proves (i). 

(ii) Using the reciprocity theorem for Gauss sums ([3]), we have for all integers a, b 
with ac ^ and ac + d even that the quantity 



S{a, b, d) := ^exp ( -jiaf + bj) 

3=0 



obeys 



S{a, b, d) 



d 






exp 1 


a 





HI 



Thus \S{a,b,d)\ = y/d if and only if \S{—d, —b,a)\ = ^Ja. Applying if to a = A; = 
\,...,d— \ coprime with d and b = 2m + k, and taking the complex conjugate yields 
the result. Namely ad + b = dk + k + 2m is even for all k = 0, ...,d — 1 since d is 
odd. 

Corollary 3.9 Let d > 3 be a prime number. Then for any k = 1,2, ...,d — 1 the 
sequences 



j=0,...,d-l 



are bi-unimodular. 
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Remark 3.10 This property is known, but has an extension in the non prime odd 
dimensions. See next section. 

It is known that the diagonahzation of VU'', k = 0, ...,d—l provides a set of (i+1 
MUB's for a prime p. Here we show that the same is true with the diagonahzation 
ofV'U, k = 0,...,d-l. 

Theorem 3.11 Let d > 2 be a prime dimension. Then the orthonormal bases 
defined by the unitary matrices 1, Pq, R, R^, ...,R'^~^ provide a set of d + 1 MUB's. 

Proof : For d = 2 this has been aheady proven in Section 2. For d > 3 (thus 
odd, since it is prime), it is enough to check that : 

(i) Pq,R'', k = 1, ...,d — 1 are unitary Hadamard matrices, together with 

(ii) Po*^^ k = l,...,d-l and {R^'fR^, l<k'<k<d-l. 

Since (i) has been aheady estabhshed, it remains to show (ii). Since i?^ = u^PoD^Pq 
we have 

P^R'' = a'^D'^P* 

which is trivially an unitary Hadamard matrix (since Pq is, a is of modulus 1 and 
D is diagonal and unitary). For (P*)^ i?^ it is trivial since 

{R*f'R^ = R^-''' 

which is unitary Hadamard for any k ^ k', k,k' = l,...,d— 1. 

4 THE CASE OF ARBITRARY ODD DIMEN- 
SION 

We have shown in [5] that for any odd dimension d > 3 the matrices P^Pk is an 
unitary Hadamard matrix provided k is co-prime with d. This can be transfered 
to a similar property for the matrix R~'^, and therefore to the bi-unimodularity of 
the sequence g^''^ defined in f ll.Sp . 

Proposition 4.1 Let d > 3 be an odd integer, and k be any number coprime with 
d. Then 

(i) R^ is an unitary Hadamard matrix. 

(ii) The sequence g^^'' is bi-unimodular. 
(Hi) Both properies are equivalent. 

This implies Proposition 1.10. 

Theorem 4.2 Let d be odd and k > 2 be the smallest divisor of d. Then the 
orthonormal bases defined by the unitary matrices {l, F, i?, i?^, i?'^"^} provide a 
set of k + 1 MUB's in dimension d. 
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5 THE CASE OF ARBITRARY EVEN DIMEN- 
SION 

Let d be even and 

/2i7r 
a; = exp ^— 

We denote cu^/^ = e*'^/'^. One defines the Discrete Fourier Transform F as usually. 
The theory of circulant matrices is also pertinent for even dimensions d > 4. Namely 
in that case the matrix 

D' = dmg{l,uj-'/\...,iu-''"/^...,u;-'/^) 

has been shown ([5]) to be such that the unitary Hadamard matrix 

Pi = D'F 

diagonalizes VU. But the circulant matrix R obtained by multiplying the k-th 
column vector of Pi by uj~'' also diagonalizes VU : 



Proposition 5.1 The circulant matrix whose matrix elements are 

1 

Vd 

diagonalizes VU and is such that F*R is an unitary Hadamard matrix. 
Proof : 

thus 

1 (i-fc)^ 

Rj^k only depends on j—k thus is circulant (and unitary). Therefore it is diagonalized 
by F, namely there exists an unitary diagonal matrix D" such that 

F*RF = D" 

Again this implies that F*R is an unitary Hadamard matrix. 

n 

Corollary 5.2 (i) The orthonormal bases defined by the unitary matrices {1,F,R} 

provide a set of 3 MUB's in arbitrary even dimension d. 

(a) One has the following property of quadratic Gauss sums for d even : 



E'^ ^ ( ik'^TC 
k=0 ^ 



Vd 
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Remark 5.3 is circulant, unitary, hut not Hadamard. Thus it does not help to 
find more than 3 MUB's in even dimensions. In dimensions d = 2", another method 
is necessary to prove that there exists d + 1 MUB's. 
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